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We study the application of decoupling techniques to the case of a damped vibrational mode of 
a chain of trapped ions, which can be used as a quantum bus in linear ion trap quantum computers. 
We show that vibrational heating could be efficiently suppressed using appropriate "parity kicks" . 
We also show that vibrational decoherence can be suppressed by this decoupling procedure, even 
though this is generally more difficult because the rate at which the parity kicks have to applied 
increases with the effective bath temperature. 



I. INTRODUCTION 



Real world quantum systems interact with their environment to a greater or lesser extent. No matter how weak the 
coupling with such an environment, the evolution of an open quantum system is eventually affected by nonunitary 
features like decoherence, dissipation, and heating. Decoherence, in particular, is a serious obstacle to all applications 
exploiting quantum coherence, such as the bourgeoning field of quantum information processing. 

Recently, considerable effort has been devoted to designing strategies able to counteract the undesired effects of the 
coupling with an external environment. Notable examples of these strategies in the field of quantum information are 
quantum error correction codes |l| and error avoiding codes [Q] , both based on encoding the state to be protected into 
carefully selected subspaces of the joint Hilbert space of the system and a number of ancillary systems. The main 
difference between the two encoding strategies is that error avoiding codes (also called decoherence-free subspaces) 
provide a passive strategy relying on the occurrence of specific symmetries in the interaction with the environment, 
which guarantees the existence of state space regions inaccessible to noise. Quantum error correction is instead an 
active strategy in which the encoding is performed in such a way that the various errors are mapped onto orthogonal 
subspaces so that they can be diagnosed and reversed. 

A simple example of decoherence-free subspace has been recently demonstrated with two trapped ions while 
error correction codes for single qubit errors has been demonstrated only in NMR quantum information processors [Q . 
The main limitation for the efficient implementation of these encoding strategies for combatting decoherence is the 
large amount of extra space resources required Q. Correcting all the possible one-qubit errors requires at least five 
qubits [^| and if fault tolerant error correction is also considered, the number of ancillary qubits rapidly increases. For 
this reason, other alternative approaches which do not require any ancillary resources have been pursued, and which 
may be divided into two main categories: closed-loop (quantum feedback) @,||], and open-loop p|-|l4|] decoherence 
control strategies. In closed loop techniques, the system to be protected is subject to appropriate measurements and 
the classical information obtained from this measurement is used for real-time correction of the system dynamics. This 
technique shares therefore some similarities with quantum error correction, which also checks which error has taken 
place and eventually corrects it. However, the main limiting aspect of feedback schemes is the need of a measurement, 
which is always inevitably subject to the limitations due to non-unit detection efficiency. In fact, only under specific 
cases (see S) it is possible to automatically correct the error without a measurement, as in quantum error correction 
codes. In open loop control strategies instead, the system is subject to external, suitably tailored, time-dependent 
drivings which are independent of the system dynamics and do not require any measurement, but only a limited, 
a priori, knowledge of the system-environment dynamics. These external control Hamiltonians are chosen in order 
to realize an effective dynamical decoupling of the system from the environment. In this way, any undesired effect 
of the environment, such as dissipation, decoherence, heating, can be eliminated in principle. The essential physical 
idea behind these open loop schemes comes from refocusing techniques in NMR spectroscopy, now routinely used to 
eliminate unwanted interactions |15|j . Nonetheless, these decoupling methods have recently attracted a large interest 
and they have been applied in many different situations, such as the inhibition of the decay of an unstable atomic 
state jL3j], or the suppression of magnetic state decoherence The general applicability of decoupling methods 

has been discussed in fl^| , while the possibility to combine decoupling techniques together with weak-strength and 
slow-switching controls has been analysed in pq] , where the conditions under which noise-tolerant, universal quantum 
control of a system can be performed with no extra space resources, have been determined. The general algebraic 
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structure behind decoupling strategies has been also analysed in Ref. |17]], where it is shown how decoupling can be 
also considered as a dynamical symmetrization with respect to a group. This more general algebraic framework has 



also provided a unifying picture for coding and decoupling noise control strategies [ [18|jl9| . In fact, when decoupling 
open loop controls are combined together with encoding into larger Hilbert spaces, fault-tolerant universal control 
of quantum systems becomes possible even with limited control resources. For example, it has been shown that 
the Hciscnberg exchange interaction is sufficient to perform universal quantum computation if appropriately encoded 
qubits are used [ pp| ; these encoded subspaces may actually be made decoherence-free if appropriate decoupling controls 
are applied in parallel fl9| , |2l"|| . 

The main drawback of open loop decoupling procedures is that the timing constraints are particularly stringent. In 
fact, the decoupling interactions has to be turned on and off at extremely short time scales, even faster than typical 
environmental timescale {full- strength/ fast- switching or quantum bang-bang controls []12|]). In fact, perfect decoupling 
from the environment is obtained only in the infinitely fast control limit (see Section 11) and it is therefore important 
to establish in a quantitative way how effective these decoupling schemes are in a realistic situation with control pulses 
with finite strength and time duration. A detailed analysis of decoupling timescales has been performed only in [|| for 
the case a single qubit in the presence of a purely dephasing environment, and in in the case of a linearly damped 
vibrational degrees of freedom. In this latter case, Ref. pT| proved that perfect decoupling can be achieved using 
extremely fast "parity kicks" , and that significant suppression of dissipation and decoherence due to the coupling 
with a zero-temperature bath is obtained as soon as the frequency of parity kicks becomes larger than the frequency 
cutoff of the environment. In the present paper we shall reconsider the model of Ref. |[lj and extend the analysis 
to the case of a finite temperature environment. The motivation for this study is twofold. First of all it will allow 
us to establish if and how thermal effects influence the decoupling strategy, that is, if temperature introduces a new 
timescale which, together with the environmental frequency cutoff, determines the effectiveness of the parity kick 
decoupling strategy. Secondly, the damped harmonic oscillator in a finite temperature bath studied in this paper well 
describes a collective vibrational mode of a chain of trapped ions, which is used as a quantum bus in linear ion trap 
quantum computers (2^] . One of the main experimental problems for quantum information processing with linear ion 
traps is just heating of these vibrational modes |p3|l , and it is therefore extremely important to establish if the parity 
kick decoupling method of Ref. jll| is able to suppress heating and decoherence in this case. 

The paper is organized as follows: In Section II decoupling strategies in general, and the parity kick method of 
Ref. [ pT| as a particular example, are presented. In Section III the dynamics of the vibrational mode in the presence 
of a nonzero temperature bath and parity kicks is analyzed in detail and in Section IV the numerical results for both 
heating and decoherence rates are presented. Section V is for concluding remarks. 

II. DYNAMICAL DECOUPLING VIA PARITY KICKS 

The starting point of decoupling techniques is the observation that even though one does not have access to the 
large number of uncontrollable degrees of freedom of the environment, it is still possible to interfere with its dynamics 
by inducing motions into the system, which are at least as fast as the environment dynamics. This indirect influence 
of the environment is obtained through the application of suitable time-dependent perturbations acting on the system 
variables only. Let us now review the main points of the decoupling technique following the lines of Ref. ]L2] ] . 

We consider a quantum system S coupled to an arbitrary bath B, whose overall Hamiltonian can be written as 

H = H s ® 1 B + Is <8> H B + H SB = ^ S a <E> B a . (1) 

OL 

A decoupling strategy consists in trying to protect the evolution of S against the effect of the interaction H$b, by 
seeking a perturbation Hi (t) <g> 1b to be added to Ho so that the total Hamiltonian becomes H (t) = Ho + -Hi (t) ® 1b- 
One usually restricts to situations where the control field is cyclic, i.e., associated to a decoupling operator U\{t) that 
is periodic over some cycle time T c : 

Ui(t) = Texp|- (i/b) J duffi(«)| = Ux(t + T c ) , (2) 

where T denotes time ordering. In this case one focus on the stroboscopic evolution at times T/v = NT C , and it is 
possible to see that in this case the evolution is driven by an effective average Hamiltonian 

U tot (T N ) = e-WVWs . (3) 
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The calculation of the average Hamiltonian H is performed on the basis of a standard Magnus expansion of the 
time-ordered exponential defining the cycle propagator f25|] , 

U tot (T c ) = exp(-iHT c /h) = Texp j- duH(u)^ = e -i[H m +H W +~ ] t c /» > (4) 

where 

ff(t) = UI^HaUxit) = [ul(t)S a Ui(tj\ ® . (5) 

The various contributions in the right hand side of Eq. (|^) collect terms of equal order in H(t), In particular, 

i? (0) = l f TC duH(u), (6) 

= ~ w r dv f du [ k{vi k[u) \ ■ (7) 

c J o J o 

One says that fcth-order decoupling is achieved if the control field Hi (t) can be devised so that contributions mixing 

S and B degrees of freedom are no longer present in and the first nonvanishing correction arises from H , 

k > 1. One then considers the infinitely fast control limit, which, for a finite evolution time T, requires considering 
T c = T/N in the limit T c — > and N — > oo. In this limit, first-order decoupling is sufficient, contributions higher 

than zeroth-order are negligible in ([I]), and one can focus on the problem of designing the effective Hamiltonian 
of Eq. (||) in such a way that there is no residual system-environment coupling. 

The more general way to engineer the average Hamiltonian is through symmetrization with respect to a finite 
group Q ]T^ , |T7| . In fact, if we consider a finite group of unitary operators Q = {gj},j = 1, • • • , \G\, symmetrization is 
the map (acting on system operator only) 

S -IMS) =1^1 9)Sg it (8) 

which is also theprojection on the so-called centralizer of Q, composed of operators commuting with every element gj 
of the group g MM. The map (^) can be dynamically implemented through a simple piecewise constant decoupling 
operator: 

U x {t)= gj , jAt<t<{j + l)At, (9) 
corresponding to a partition of the cycle time T c into \Q\ intervals of equal length At = T C /\Q\. Then, by (||), 

H {0) - Uc(H ) = u c(S a ) ® B a , (10) 

a 

showing that the average Hamiltonian , generating time evolution in the infinitely fast control limit, has been 
symmetrized, i.e., has become invariant with respect to the group Q. Perfect decoupling from the environment is 
achieved when 

Tl c (H SB ) = 1 (11) 

and in this case the effective open system evolution for the reduced density operator of the system ps over time T is 
governed by 

lim Ps (T = NT C ) = e - t7IsT / h p s (0) e +l ~ EsT/h , (12) 

where Hs = Tlc(Hs). This means that the system is no more interacting with the environment and the residual time 
evolution is driven by a projected system Hamiltonian, invariant with respect to Q . 

A number of examples of decoupling groups Q has now appeared in the literature, expecially for the case of many- 
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qubits dissipative registers with various kinds of interaction with the environment |]9|,[lO| Jl^Jl7| ]. Another important 
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example for applications in quantum computing is the case of a linearly dissipative vibrational degree of freedom, 
which can be used as a quantum bus in linear ion trap quantum computers [^2| , and which has been shown in Ref . Jll]] 
to be decoupled by the group Z2, composed by the identity and the parity operator P. In fact, it is straightforward 
to check that the decoupling condition (|ll|) is equivalent to the condition of Eq. (7) of Ref. jlj)] . 

As it can be expected, implementing the above general decoupling strategy is by no means trivial. First of all, 
for a given H$b, the identification of a minimal group Q able to produce decoupling is nontrivial. Secondly, the 
decoupling prescription (^) requires the capability of instantaneously changing the evolution operator from gj to gj+i 
over successive subintervals. This means assuming the capability of implementing arbitrarily strong and extremely 
fast control operations. Such impulsive full-power control configurations correspond to so-called quantum bang-bang 
controls as introduced in (|] . As it has been already shown in J9|,[ll| , the most stringent condition is not on the strength 
but rather on the extremely high speed of the control operations: one has to be faster than the typical environmental 
timescale, which is usually fixed by the frequency cutoff of the bath spectrum, lu c . However, the identification of the 
frequency cutoff to c as the only relevant parameter determining the threshold for the decoupling cycle frequency 1/T C 
above which the decoupling procedure becomes effective, has been done in Refs. only on the basis of two specific 

examples. In Ref. |n|] , the case of a harmonic oscillator coupled to a zero-temperature bath, able to induce only 
system dissipation (and the associated decoherence) has been considered. The case of nonzero temperature has been 
discussed in Ref. |j) but only in the particular case of a single qubit subject to a purely dephasing, energy-conserving, 
environment. It is therefore important to establish the effectiveness of decoupling techniques in the general case of a 
nonzero-temperature, dissipative bath. From now on we shall specialize to the case of the linearly damped harmonic 
oscillator of Ref. which is of relevance for linear ion trap quantum computation. In fact, we shall demonstrate 
that decoupling techniques can be successfully used to efficiently suppress heating of the vibrational center-of-mass 
motion of the ion chain. 



III. PARITY KICKS FOR A DAMPED HARMONIC OSCILLATOR 



We choose a harmonic oscillator as system of interest 

Hs = Titu^a) a , (13) 

describing a collective vibrational mode of a linear chain of trapped ions with frequency loq. It has been already 
experimentally verified | f23|]2(|| that the nonunitary features of the vibrational dynamics (heating and decoherence) 
are well described by modelling the environment as a collection of independent bosonic modes |^7j 

H B =J2 M{b k , (14) 

k 

interacting with the vibrational mode via the following bilinear interaction Hamiltonian in which the "counter- 
rotating" terms are dropped 

H SB =J2 n, yk(abl+a%k) ■ (15) 

k 

The symmetrization with respect to the group Z2 is performed by periodically pulsing the oscillation frequency, that 
is, by changing the potential so that ujq is changed to luq + Suj for a time interval t and with a time period T c (see 
Fig. 1). The pulse realizes the "parity kick" of Ref. JllJ] when the condition Slu ■ r = ir is satisfied. In this way, the 
cyclic time-dependent control Hamiltonian is given by 

00 

#i(i) = hSua^a^O (t - nT c + r)) 9 {nT c - t) , (16) 

n=l 

so that the cyclic decoupling operator U\{t) is equal to U\{t) = Is for < t < T c — r and U\{t) = expliira^a} = P, 
for T c - t < t < T c . 



4 



co + 5(0 



co 



FIG. 1. Sketch of the implementation of the parity kick decoupling procedure by pulsing the oscillation frequency. 



To determine the effects of a nonzero temperature bath on the efficiency of the decoupling scheme, we shall consider 
an initially factorized state in which the vibrational mode is prepared in a given pure state |^(0)) and the environment 
is at the thermal equilibrium state at temperature T, To be more specific we want to establish if decoupling 
via parity kicks is able to suppress efficiently both heating (which is important for quantum information processing) 
and quantum decoherence of the vibrational mode. To study heating we shall assume that the collective vibrational 
mode has been initially cooled to its ground state p8| , that is, \4>{0)) = 0). The study of decoherence instead will be 
performed, as in Ref. UM, by considering an initial linear superposition of two coherent states with opposite phases, 
that is, the well known Schrodinger cat state 



mo)) = Wv> = N <p (Ho)) + e lv \ - «(o)» 



(17) 



where N v = (2 + 2e~ 2 l"(°)l costp) -1 / 2 . The dynamics of the system in the presence of parity kicks will be exactly 
solved for both initial conditions of the vibrational mode, by exploiting the fact that a tensor product of coherent 
states retains its form at all times when the evolution is generated by the Hamiltonian of Eqs. (|l4|) and ( |T5| ) 

PI, that is 



«(o)> ® n i&(°)> -> i«(*)> ® n i&(*)> . 



k k 

where the time-dependent coherent state amplitudes are a linear combination of the initial amplitudes 



a(t) = Loo(*)a(0) + ^ L ak (t)p k (0) 

k 

0k(t) = L kQ (t)a(O)+J2L k k'(t)0k'(O). 



(18) 

(19) 
(20) 



Eq. (]l^) is useful also in the nonzero temperature case. In fact, using the expression of the thermal state p\ in the 
Glauber-Sudarshan P-representation pel] 



T 

Pb 



n 



d 2 (3 k 
irN k 



exp 



\Pk)(Pk 



(21) 



where d 2 (3t = dRe/3fc<iIm/3fc and N k = (exp {hujk/kBT} — 1) _1 is the mean thermal excitation number of the k-th bath 
mode, one has always to evaluate the time evolution of terms like |a)(a'||/3fc)(/3fc|, and then perform the average over 
the thermal Gaussian weight exp { — \(3k\ 2 /Nk} JjrNk- Therefore the essential dynamics is contained in the expression 
of the unitary matrix Ly(t) of Eqs. ( JToj ) and (pp[), which has the same structure both with and without parity kicks, 
because the two situations differ only by the value of the oscillation frequency. The matrix element Loo {t) is given in 
terms of its Laplace transform, and, in the interaction picture with respect to Hs of Eq. (O), one has 



where 



L 0O (t,8u) = C- 1 



if(z,M = E 



z + K(z, Suj) 



it 



z + i(u>k — luq) — iSlu 



(22) 



(23) 
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This expression refers to the evolution during the parity kicks, that is, for nT c — r < t < nT c , n > 1. The evolution 
in the absence of kicks is simply obtained putting Slu — in Eqs. ((2^) and All the other matrix elements can be 
expressed in terms of the matrix element £oo(i; 5u>) in the following way: 

L ok (t, Su) = L M (t, Su) = -*7fc / d Se - l ^ fc - Wo -^) s L o(i - s, <M (24) 



L kk ,(t,Suj) = 5 kk ,e~ l ^ k ~^ a)t - lkl ' k / d 8e -<(w»-«o-&.')(*-<0 / d s ' e - l K'-^-' 5 ")( ;i - s ')Loo(s',fe) • (25) 

■/o Jo 

It is evident that a decoupling cycle of duration T c will be described by the product of unitary matrices L(t, Slo) ■ 

L(T C — r, 0) applied to the vector formed by the coherent amplitudes (a(t), Pk(t) )• As a consequence, the 

stroboscopic dynamics of the whole system during the decoupling procedure, in the case of an initial tensor product 
of coherent states, can be described exactly as 



Pk{NT c ) 



\ 



[L(t,6u)-L(T c -t,0)] 



N 



) 



I a(0) ^ 

V ; / 



(26) 



IV. NUMERICAL RESULTS 



In the standard description of dissipation, one always considers a continuum distribution of oscillator frequencies in 
order to obtain an irreversible transfer of energy from the system of interest into the reservoir. Moreover, most often, 
also the Markovian assumption is made which means assuming an infinitely fast bath with an infinite frequency cutoff 
lj c . This case of a standard vacuum bath in the Markovian limit is characterized by an infinite, continuous and flat 
distribution of couplings ||(]] , 

7H 2 = _L Vw , (27) 

where 7 is the energy damping rate. As shown in Refs. |j 11 1, decoupling strategies become efficient when the 
external controls are characterized by timescales faster than those of the environment. It is therefore evident that, 
in the presence of parity kicks, we cannot make any Markovian approximation. We have to solve numerically the 
problem, by simulating the continuous distribution of bath oscillators with a large but finite number of oscillators 
with closely spaced frequencies. As in Ref. |ll|, we have considered a bath of 201 oscillators, with equally spaced 
frequencies, symmetrically distributed around the resonance frequency i.e. 

Lu k =uj + kA A = ^ (28) 
k max = ^ = 100 ior x = 2^o (29) 

kmin = -kmax = -100 => uf in = , (30) 

(31) 

and we have considered a constant distribution of couplings similar to that associated with the Markovian limit 

ll = ^ Vfc . (32) 

Approximating a continuous Markovian bath with a finite number of bath oscillators has two main effects. First of 
all, the discrete frequency distribution with a fixed spacing A makes all the dynamical quantities periodic with period 
T rev = 2ir/A Therefore our numerical solution will correctly describe the interaction with the environment 

provided that we consider not too large times, say t < 7r/A. Secondly, the introduction of a finite cutoff (u> c = 2wo in 
our case) implies a modification of the coupling spectrum 7(0;) at very high frequency with respect to the infinitely 
flat distribution of the Markovian treatment (see Eq. (|27])). This fact manifests itself in a slight modification of the 
dynamics at very short times (t ~ a;" 1 ) pl| . We have verified both short and long time deviations from the standard 
Markovian bath dynamics in our numerical calculations. However, we have checked that our model environment with 
a finite number of oscillators faithfully reproduces the standard Markovian bath dynamics within the time interval of 
interest, O.I/7 < t < 3/7 say. 
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A. Effect of parity kicks on heating 



To check if decoupling via parity kicks is able to suppress the heating of the vibrational mode, we consider the 
following initial state for the whole system 

i°><°i®n / ^ ex p{-^F}i^(°))^(°)i' ^ 

where |0) is the ground state of the collective vibrational mode p8| . 

Using the property (|l^), and tracing over the environment, the evolved state of the vibrational mode after N 
decoupling cycles can be written as 

Ps(NT c ) = /n%P ex p{"^f^} HNT c ))(a(NT c )\, (34) 

where the coherent state amplitude a(NT c ) is the following linear combination of the complex variables /3fc(0), 

a(NT c ) = Y, C ok(NT c )p k (0), (35) 

k 

where we have defined Cok{NT c ) = ([L(t,Slj) ■ L(T c — t,0)] N X (see Eqs. ( |i~9| ) and (p6|)). The Gaussian average 

of Eq. ( |34| ) can be performed by first considering the normally ordered characteristic function x(A,iVT c ) Q of the 
state, and then performing the integration. One gets 

x(A,7VT c ) = exp|-|A| 2 ^7V fe |c ofe (7VT c )| 2 |, (36) 

showing that, in the presence of parity kicks, the vibrational state is a thermal state, with mean vibrational number 
v{NT c ), 

v{NT c ) = J2 N k\Cok(NT c )\ 2 . (37) 

k 

The stroboscopic time evolution of this mean vibrational number is plotted in Fig. 2 both in the presence (full 
circles) and in the absence (crosses) of parity kicks. The capability of the parity kick decoupling strategy to avoid 
vibrational heating is clearly visible in this figure. In Fig. 2 and in the rest of the paper we consider a vibrational 
mode with frequency loq = 10 Mhz, damping rate 7 = 0.1 Mhz and environmental frequency cutoff lo c = 20 Mhz. 
The curve referring to the situation without parity kicks in Fig. 2 well reproduces the standard Markovian result |3(]] 
v{t) = N(u)q)(1 — e~ 7 *), where N(uj ) — (exp {ftwo/fceT} — 1) _1 is the mean vibrational number of the oscillator 
at thermal equilibrium in the usual Born-Markov approximation. Fig. 2 refers to an effective reservoir temperature 
T = 10 mK (corresponding to N(ui ) ~ 130), and to the following decoupling cycle parameters: T c = 157 ns, parity 
kick duration r = T c /7 ~ 22.4 ns, implying Suj = 140 Mhz. 
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FIG. 2. Time evolution of the mean vibrational number of Eq. (|37j) with (full circles) and without (crosses) parity kicks. 
The capability of parity kicks to suppress heating is clearly visible. Parameters are: uo — 10 Mhz, 7 = 0.1 Mhz, u> c = 20 Mhz, 
effective reservoir temperature T = 10 mK (corresponding to N(u)o) — 130), T c = 157 ns, parity kick duration r = T c /7 ~ 22.4 
ns, implying Suj — 140 Mhz. 



The influence of the environmental temperature T on heating suppression is analysed in Fig. 3, where the mean 
vibrational number after one relaxation time t — 1/7, 2/(1/7), IS plotted as a function of the rescaled decoupling 
cycle time uj c T c /2n for three different bath temperatures, T = 10 mK (a), T = 100 mK (b), and T = 1 K (c). For 
each value of T c , we have always chosen the kick duration r = T c /7, as in Fig. 2, and the frequency shift Sco is 
always correspondingly adjusted so that Sui = tt/t. We can see that a well visible threshold for the decoupling cycle 
time T c exists and that as soon as the parity kicks are sufficiently fast, T c < 2t:/lu c , heating suppression becomes 
significant. This is a p hase transition- like behavior analogous to that found for decoherence suppression in the zero- 
temperature case [[Ll|. What is more important is that bath temperature has no effect on the effectiveness of the 
decoupling scheme: the results are essentially identical for the three different temperatures studied, and this means 
that, at least for what concerns heating, the only relevant environmental timescale is given by the frequency cutoff 
lu c . This result is particularly important for the application of the parity kick strategy to suppress heating in linear 
ion trap quantum computers, where heating is due to some technical imperfections originating from fluctuating patch 
fields p| . Determining the effective temperature T of the thermal bath modelling these fluctuating fields is generally 
very difficult, but our results shows that this is not relevant, and that parity kick decoupling is very promising for 
eliminating vibrational heating. 
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FIG. 3. Mean vibrational number after one relaxation time t = I/7, 2/(1/7), as a function of the rescaled decoupling cycle time 
lu c Tc/2h, for three different bath temperatures: T = 10 mK (corresponding to N(loq) ~ 130) (a), T = 100 mK (corresponding 
to N(loo) ~ 1302) (b), and T — 1 K (corresponding to N(ujo) — 13144) (c). For each value of T c , we have always chosen 
r = T c /7, and, correspondingly, Scu = tt/t. The other parameters are as in Fig. 2 
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B. Effect of parity kicks on decoherence 



Let us now consider the possibility to suppress decoherence. We assume an initially prepared Schrodinger cat state 
of the vibrational mode, and therefore the following initial state for the whole system 

hM<vvi®n / %r ex p{^^l^} l/3fc(0))(/3fc(0)l ' (38) 

where W) 1S given by Eq. (|l7|). 

Using the property (|l^), and tracing over the environment, the evolved state of the vibrational mode after N 
decoupling cycles can be written as 

p s (NT c )=N* J n^^exp|-^M|{|a + (7VT c ))(a + (iVr c )| + |a_(7VT c ))( a _(iVr c )| 

W*(/?+(iVT c )|/3 fe -(Ar^ (39) 

The coherent state amplitudes a±(NT c ) and f3^(NT c ) are now given by the following linear combinations of the 
initial amplitudes (see Eqs. (^9|)-(p0|)) 

a±(NT c ) = ±a Q C QQ (NT c ) + ^ C ok {NT c )(3 k {Q), (40) 

k 

pf(NT c ) = ±a C k o(NT c ) + £ C kk ,(NT c )0 k ,(O), (41) 



The Gaussian average of Eq. (|39|) can be performed, as in the preceding subsection, by first considering the normally 
ordered characteristic function of the state and then performing the integration. The integration is straightforward 
but lengthy, and the resulting reduced vibrational state may be better expressed in terms of its Wigner function 

W s (a,NT c ), 



gjVg f f 2\a-a a C oa (NT c )\ 2 \ f 2\a + a C 00 (NT c )\< 

7T [1 + 2u{NT c )\ \ CXP \ 1 + 2u{NT c ) / exp \ 1 + 2v{NT c ) 

4Im[aaoCoo(iVTc)]" 



+2exp{-2|a^(iVT c )}exp{- IT ^L} 



cos 



1 + 2v{NT c ) 



(42) 



where v(NT c ) is again the mean vibrational number of the cat state of Eq . (p7|), the matrix element Cqo(NT c ) 
describes the amplitude decay, and r}{NT c ) is the fringe visibility function |32|, determining the relative strength of 
the quantum interference term in the cat state, and which can be expressed as 

n (NT c ) = l- \ C ^ NT o)? (43) 

This fringe visibility is always contained in the interval [0, 1] and provides a good quantitative description of dynamical 
decoherence processes. For this reason we shall study the stroboscopic evolution of this quantity, as in Ref. pT| , to 
quantify the eventual decoherence suppression caused by the decoupling. 

The time evolution of the fringe visibility is plotted in Fig. 4, both with (full circles) and without (crosses) parity 
kicks. The possibility to suppress decoherence using parity kicks is clearly demonstrated in this figure. Parameters 
are the same as in Fig. 2, except that the decoupling cycle parameters now are: T c = 78.5 ns, parity kick duration 
t = T c /7 ~ 11.2 ns, implying 8lj = 280 Mhz. The curve referring to the situation without parity kicks in Fig. 4 
(crosses) well reproduces the Markovian result Q 

V(t) = 1 — — — — (44) 

/w l + 2iV(wo)(l-e-T*) V ' 
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FIG. 4. Time evolution of the fringe visibility function of Eq. (]43|) with (full circles) and without (crosses) parity kicks. The 
capability of parity kicks to suppress decoherence is clearly visible. Parameters are the same as in Fig. 2, except that the 
decoupling cycle parameters now are: T c = 78.5 ns, parity kick duration r = T c /7 ~ 11.2 ns, implying 8uo — 280 Mhz. 

The influence of temperature on decoherence suppression is studied in Fig. 5, where the fringe visibility function 
after one relaxation time t — 1/7, 77(1/7), is plotted as a function of the rescaled decoupling cycle time u> c T c /2n for 
three different bath temperatures, T = 10 mK (a), T = 100 mK (b), and T = 1 K (c). For each value of T c , we have 
always chosen the kick duration r = T c /7, as in Figs. 2 and 3, and the frequency shift 8lj is always correspondingly 
adjusted so that Su> = tt/t. 

We can see from Fig. 5 that the situation is rather different from that with heating suppression. In fact, decoherence 
suppression by parity kicks strongly depends on the bath temperature, and it is significant only in the lower temperature 
case (Fig. 5a), which is the only case in which a threshold for the decoupling cycle time T c at about T c ~ 2ir/u> c , as 
in the zero temperature case fllf , is visible. In the other cases, decoherence suppression worsens for increasing bath 
temperature. This result shows that eliminating decoherence via decoupling techniques is generally more difficult 
than eliminating heating. This can be easily explained in terms of the so-called thermal acceleration of decoherence 
p2|,p3[ , that is, the fact that in the case of a thermal bath at temperature T, the decoherence process is accelerated 
roughly by a factor (1 + 2iV(co>o)) with respect to the zero temperature case. This thermal effect on the decoherence 
rate can be also easily checked from the Markovian limit expression of Eq. (0). In fact, the fringe visibility function 
T](t) reaches its asymptotic value in a time of the order of td ec — [7(1 + 2-/V(wo))] -1 , and it is evident that decoherence 
suppression with parity kicks is possible only if the cycle time T c is smaller than this decoherence time tdec, and not 
only smaller than 2tt/uj c , as in the zero temperature case. This means that, in a nonzero temperature bath, one has 
a new, temperature- dependent, threshold for decoherence suppression, given by 



T c < min 



|2^K,[7(l + 2iVM)]- 1 }, 



(45) 



and this generalized expression easily explains the results of Fig. 5. 
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FIG. 5. Fringe visibility function after one relaxation time t = I/7, 77(1/7), as a function of the rescaled decoupling cycle time 
uj c T c /2ti, for three different bath temperatures: T = 10 mK (corresponding to N(luo) ~ 130) (a), T = 100 mK (corresponding 
to N(ujo) ~ 1302) (b), and T = 1 K (corresponding to 7V(wo) — 13144) (c). For each value of T c , we have always chosen 
t = T c /7, and, correspondingly, <5cj = tt/t. The other parameters are as in Fig. 2. The quality of decoherence suppression 
degrades with increasing temperature. 



V. CONCLUSIONS 



We have studied the application of open-loop decoupling schemes in an experimentally realistic scenario. In fact, 
decoupling strategies have been proved to provide perfect isolation of a system from its environment in the infinitely 
fast control limit, i.e., in the case of very intense and very fast control pulses jl2). The efficiency of decoupling 
strategies in concrete situations involving finite strength and finite duration control pulses has been analysed only in 
the specific cases of a sing le qubit in a nondissipative environment in Q , and for a damped harmonic oscillator in a zero- 
temperature bath in jll| . Here we have extended these studies to the case of a dissipative and nonzero-temperature 
reservoir. We have specialized to the case of a collective vibrational mode of a linear ion chain, which is used as a 
quantum bus in linear ion trap quantum computers |22j |. We have shown that the parity kick decoupling strategy 
introduced in JlT[ can be successfully applied to su ppre ss vibrational heating, which is one important limitation for 
quantum information processing in linear ion traps p3fl . In fact heating is suppressed as soon as the decoupling cycle 
time T c becomes smaller than 27r/w c , where to c is the bath frequency cutoff, and more importantly, the efficiency 
of this suppression is not affected by the temperature of the bath. The parity kick method can be applied using 
present technologies and its experimental implementation in the case of trapped ions would be the first example 
of the application of decoupling techniques outside the field of NMR, where the so-called "refocusing" techniques 
fl5f are easier to use because the involved magnetic environment is usually very slow (see however Ref. |34|| for a 
proof-of-principle demonstration of quantum bang-bang control in a photon polarization qubit). 

We have also shown that, differently from heating, the suppression of vibrational decoherence is more difficult, 
because in a nonzero temperature bath, the threshold for the decoupling cycle frequency is determined not only by 
the bath frequency cutoff, but also by the decoherence rate, which increases for increasing temperatures. The parity 
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kick cycle frequency has to be larger than both rates and this makes suppression of vibrational decoherence more 
difficult for higher temperatures. 
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